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SUMMARY 

Microbond tests composed of single fibre and matrix droplet are often used to deter- 
mine the properties of fibre reinforced composites. Interfacial shear strength is quantified 
by the maximum pull-out force assuming a uniform stress distribution along the fibre. 
Here, nonlinear finite element analyses are performed to investigate the validity of this 
assumption. 
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INTRODUCTION 



When determining the interfacial shear strength by the microbond test, it is commonly 
assumed that the relationship between pullout force and embeddment length is linear. 



However, many experimental results are in disagreement with this assumption [14|, [8J], 
which might be explained by the influence of boundary conditions, by the presence of 
thermal stresses and by the fracture propagation along the fibre. In this work, these 
possible effects are investigated for microbond systems consisting of a Polypropylene (PP) 
matrix and a glass fibre. 

The behaviour of fibre reinforced composites is strongly influenced by the interface prop- 
erties between fibre and matrix. One of the commonly used methods to determine these 
properties is the microbond test. In this test, a fibre embedded in a droplet of matrix 
material is pulled through a gap formed by two knife blades, so that the droplet is sheared 
off the fibre [To|, 0] - 

From this test, the average shear strength (ITTZ) is determined as 



the maximum force that is carried by the fibre divided by the area of embeddment of the 
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fibre. This average shear strength is then used for the modelling of the response of fibre 
reinforced composites consisting of many fibres. 

Although the concept of the microbond test is simple, its loading setup results in a complex 
nonuniform shear stress distribution along the fibre, which changes during the pull-out 
process. Phenomena influencing this shear stress distribution are the presence of radial 
stresses and the fracture process during pullout. Radial stresses are caused by the knife 
blades pressing against the droplet and the nonuniform geometry of the droplet itself. 
In addition, radial stresses are also generated by the difference in thermal contraction 
of fibre and matrix during the manufacturing process of the test setup, which involves 
cooling down of the fibre and the melted matrix. 

In this study, the microbond test is analysed by a three-dimensional finite element model. 
The fracture process is described by nonlinear fracture mechanics (NLFM). NLFM de- 
scribes the nonlinear response within the fracture process zone by means of softening 
stress-strain relationships. For local constitutive models, the stress at a point depends on 
the history of this point only. Fracture is represented by localised strain profiles within 
finite elements, whereby the integral of the energy dissipation density over the size of 
the finite element results in the dissipated energy. Thus, for local softening constitutive 
models for fracture, the dissipated energy depends on the size of the finite element. How- 
ever, the energy dissipated during the fracture process is a material property and should 
be independent of the discretisation applied. This limitation of local constitutive models 
can be overcome by adjustment of the softening modulus of the softening branch of the 
stress-strain curve with respect to the finite element size , Bi| . With this approach, 
the dissipated energy is modelled mesh independently, as long as the inelastic strains lo- 
calise in a zone of assumed size. This approach, which was used earlier for the analysis 
of delamination in sandwich structures [if and cracking in cohesive materials Q, is used 
in the present study. 

The present NLFM approach differs strongly from linear elastic fracture mechanics (LEFM) 
approaches, which assume a large stress-free crack with all the nonlinearities of the frac- 
ture process concentrated in an infinitesimally small zone in front of the crack tip 
If the size of the fracture process zone is large compared to the length of the stress free 
crack and the size of the structure, LEFM results in a poor approximation of the fracture 
process. It is believed that this is the case for the polyprolyene matrix system. There- 
fore, nonlinear fracture mechanics (NLFM) is chosen, which provides a more accurate 
description of the fracture process in these situations. 

The aim of this work is to investigate the influence of radial stresses and the fracture 
process on the average shear stress obtained from the microbond test. A parameter study 
is carried out to quantify, which of the influences is dominant for the case of the microbond 
test. To the authors' knowledge, nonlinear fracture mechanics has not been used before 
to analyse the microbond test. It is hoped that with this study, further insight into the 
microbond test is gained. 
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FINITE ELEMENT MODEL 



The present section deals with the description of the geometry, boundary conditions and 
material laws used in the finite element modelling of the microbond test. 



Droplet geometry and boundary conditions 



The shape of the droplet is governed by the way it is manufactured. This involves the 
heating of a small amount of matrix material, attached to the fibre, until a droplet of 
the melted matrix material forms around the fibre. The shape of the droplet is mainly 
determined by attractive forces between matrix and fibre, since the volume of the droplet 
is so small that gravitational effects are negligible. The droplet geometry is derived by the 
Laplace excess pressure across the droplet surface as discussed in 0, ^ . The geometry of 
the fibre droplet composite (Fig. [T^, and b) is determined by the fibre radius rj, the droplet 
diameter d and the embedded length t. In the finite element model shown in Fig. [Tb, only 
a quarter of the droplet is modeled. Furthermore, an interfacial transition zone (ITZ) of 
thickness h e is introduced. This three-dimensional representation allows one to consider 
the influence of the knife blade. This representation is advantageous over axisymmetric 
finite element models, which automatically assume axisymmetric boundary conditions. 
In the literature, results of modified microbond tests have been reported, where the fibre 



is pulled through a circular hole in a plate [141 ] . This modified setup results in boundary 
conditions which are closer to axisymmetric. 

The fibre end is subjected to a force F in x-direction. The constraint of the droplet in 
the x-direction due to the knife blade is modeled by a kinematic boundary condition in 



3 



Table 1: Model parameters of the structural analyses. All the values are given in /im. 

Droplet size I d r h e e c 
Small 100 27 8 0.5 10 
Medium 215 77 8 0.5 10 
Large 400 164 8 0.5 10 



one node (Fig. [Tb). The gap between the knife blade and the fibre is denoted as e c . 
In the nonlinear finite element analysis, the pull-out force was controlled by an indirect 
displacement control using the relative displacement between the droplet and the fibre at 
the boundaries of the interfacial transition zone. 

One of the key points of this study is the influence of the embedded length on the interface 
strength and how this parameter is affected by radial stresses and the fracture process 
zone. To be able to investigate this effect, three different droplet sizes were analysed with 
their dimensions presented in Tab. [TJ In the analyses for the three specimen sizes, the 
fibre radius rf and the knife gap e c were kept constant for all three specimen sizes (Tab. [I]). 

Constitutive model 

The material response of the fibre and the droplet is described by small strain isotropic 
linear elasticity with the stress-strain relationship 

<r = D:e (1) 

where cr is the stress tensor, e is the strain tensor and D is the isotropic linear elastic 
stiffness tensor based on the Young's modulus E and Poisson's ratio v. 

The nonlinear interaction between the droplet and the fibre in the ITZ is modelled by 
an isotropic elasto-damage model in a thin layer of continuum elements in the droplet 
adjacent to the fibre. The isotropic continuum model has the stress-strain relationship 

cr = (I - uj)T) : £ = (I - uj) cr (2) 

Here, uj is the damage variable and cr is the effective stress tensor, i.e. the stress that 
is carried by the undamaged material. The damage variable uj is a function of history 
variable k which is defined as the maximum equivalent strain 7 reached in the history of 
the material: n(t) = max7(r) for r <t, where t is the time. 

The equivalent strain is defined so that the Mohr-Coulomb strength envelope (Fig. [2^) of 
the form 

|r| + <7tan0 — c = (3) 

is obtained for 7 = 70 = (ccos0)/G, where c is the cohesion, i.e the shear resistance at 
zero normal stress a, G is the shear modulus and <fr is the angle of internal friction. 
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(a) 
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Figure 2: a) The initial stress envelope represented in the principal stress space obtained 
with the damage loading function, b) The softening stress-strain curve obtained with the 
damage model. 



According to Fig. [2^ and 0, chap. 15.1.5], the strength envelope in Eq. ([3]) can be written 

c ^— — -tan0Jcos0 (4) 



as 

01 — <7 3 / 0-i+ 3 



This gives, 

(1 + sin 0) 0i — (1 — sin 0) 3 = 2c cos (5) 
For a pure shear state, r m = oi = — 3 and Eq. (jSJ) reduces to 

2r m = 2c cos (6) 

The equivalent strain in this pure shear stress state is chosen to be 7 = r m / G. For general 
stress states, the equivalent strain definition is determined by setting Eq. (jSJ) and Eq. (jSJ 
equal and solving for 7. This gives, 

7= 2^((1 + sin 0)cri - (1- sin 0)0-3) (7) 



The strength envelope in Eq. turns into the Tresca strength envelope for = 0° and 
into the Rankine strength envelope for = 90°. For the Tresca strength envelope, a 
pure shear stress state is limited by the cohesion c independent of normal stresses. For 
the Rankine criterion, for which the shear strength is strongly dependent on the normal 
stress, the strength is determined by the maximum principal stress component. These two 
limiting cases will be used to investigate the influence of thermal stresses on the response 
of the microbond test. For = 0°, the interfacial strength given by the microbond test is 
expected to be independent of the thermal stress. However, for = 90° a strong influence 
is expected. 
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Table 2: Material parameters for the three phases at a temperature of T = 20 °C. 



Phase 


E [GPa] 


V 


c [MPa] 


G c [J/m 2 ] 


a T [1/K] 


Fibre 


7.6 


0.22 






1 x 10~ 6 


Matrix 


2.5 


0.2 






111 x lO" 6 


ITZ 


2.5 


0.2 


5 


3000 


1 x 10~ 6 



The damage variable to is related to the history variable k as 

f if«<7o 

u = g («) = < 7o / K-7o \ ., . (8) 
1 exp if k > 70 

I « V 7f-7o/ 
The parameter 7f is related to the fracture energy G*f in shear as 

(jrp 1 . . 

ll= E^h c + 2 7 ° (9) 

where /i e is the depth of the element row (softening band) next to the fibre, which rep- 
resents the ITZ. This damage law results in a softening stress-strain curve in shear, as 
presented in Fig. Eb. The fracture energy, which is expected to influence strongly the 
interfacial shear strength, is one of the parameters, which will be varied in the present 
study. The material parameters of the three phases were chosen according to Tab. [2] for 
a temperature of T = 20 °C. 

The matrix material undergoes a phase change during the cooling down, which influences 
the mechanical and thermal properties. This is considered in the constitutive model by 
temperature dependent functions for the Young's modulus and the thermal expansion 
coefficient. A nonlinear thermal expansion coefficient of the form 

a T (T) = a t T + b t (10) 

was chosen. Furthermore, the temperature dependent stiffness of the matrix was chosen 

as 

E (T) = E exp(-((T + a c )/b e ) 2 ) (11) 
The parameters in Eqs. (llOp and (II ip are determined by a best fit to experimental results 



reported in [13] as a t = 9.9 x 10" 7 , b t = 9.2 x 10~ 5 , £ = 5 GPa, a e = 75 °C and b e = 75 °C 
(Fig. [3]). For the fibre and ITZ, the thermal expansion coefficient and Young's modulus 
in Tab. [2] for T = 20 °C was assumed to be independent of the temperature. 



RESULTS 



The results are divided into two parts. Firstly, the effect of radial stresses are investigated. 
Secondly, the influence of the fracture process zone is studied. 
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Figure 3: Comparison of thermal properties of the matrix used in the constitutive model 
to experimental results reported in [ITJ. Influence of the temperature on (a) thermal 
expansion and (b) Young's modulus. 



Cooling down of the droplet is expected to generate large radial stresses. Furthermore, 
the pressure of the knife blades and the shape of the droplet will cause radial stresses 
as well. To be able to study the effect of thermal stresses, knife blades and geometry 
separately, analyses with cooling down and without cooling down are performed for the 
three droplet sizes presented in the previous section. For each of the six analyses, the 
average interfacial shear strength is determined by dividing the maximum pullout force 
by the embedded area of the fibre. The Rankine criterion is used for the analyses, which 
results in the greatest increase of shear stresses due to radial compressive stresses. The 
results of the analyses are presented in Fig. H] in the form of the average shear stress 
normalised by the local shear strength c (Tab. [2]). 

For the thermo-mechanical analyses the average shear strengths are up to 25 % greater 
than those obtained for the mechanical analyses alone. This indicates that the cooling 
down activates significant radial compressive stresses. To investigate further the influence 
of thermally induced radial stresses and the role of the frictional angel 0, a parameter study 
was performed, in which thermo-mechanical analyses were carried out for the smallest 
specimen (100 /xm) with the frictional angle varying from 0° to 90°. The results in the 
form of the normalised average shear strength versus the frictional angle is presented in 
Fig. El These results demonstrate that the frictional angle has a strong influence on the 
increase of shear stresses due to the radial compressive stresses. In experiments reported 
in 2(|, the frictional angle for polypropylene is reported to be in the range of 10° to 
15°. For an angle of friction of 15°, the present analysis predicts the shear stress to be 
6 % greater than for = 0°. The influence of the cooling down is further investigated 
by studying the stress distribution in the ITZ along the fibre. For a frictional angle of 
cf) = 90° and an embedded length of I = 100 /j,m, the radial and shear stresses (see Fig. [Tb 
for a definition of the cylindrical coordinate system) are presented in Fig. [U] after the 
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Figure 4: Normalised average shear stress versus embedded length for thermo-mechanical 
and mechanical analysis with the Rankine criterion (0 = 90°). 
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Figure 5: Normalised average shear strength versus frictional angle <fi for thermo- 
mechanical analyses of the smallest specimen with an embedded length of t = 100 /im. 
The average shear strengths are normalised by a local shear strength of 5 MPa. 
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Figure 6: Contour plot of (a) radial stresses and (b) shear stresses in the ITZ after cooling 
down for a frictional angle of = 90° and a embedded length of £ = 100 /im. 
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Figure 7: Contour plot of (a) radial stresses and (b) shear stresses in the ITZ at maximum 
pullout force for a frictional angle of <fi = 90° and a embedded length of £ = 100 /im. 



cooling down and in Fig. [7] at the maximum pullout force. 

After the cooling down (Fig. [6]), radial compressive stresses of up to 50 MPa are present. 
However, these stresses are reduced during the shearing off, so that the maximum radial 
compressive stress at maximum pullout force is only 7 MPa. 

In the second part of the analyses, the influence of the fracture energy is investigated. 
The responses of the three droplet sizes with £ = 100, 215 and 400 /im are investigated 
for fracture energies G{ = 30, 300 and 3000 J/m 2 . The pullout process is analysed for 
a constant temperature of T = 20 °C without taking into account the cooling down 
from T = 120 to 20 °C. In this way, only the mechanical response is considered, so that 
the influence of the fracture energy can be separated from the effect of thermal stresses. 
The results are presented in the form of the normalised average shear stress versus the 
embedded length in Fig. [HJ With decreasing fracture energy, the averaged shear strength 
decreases. Furthermore, the average shear strength decreases with increasing embedded 
length. The influence of the embedded length is most pronounced for a fracture energy of 
30 J/m 2 . However, this fracture energy is much less than fracture energies of PP reported 



in the literature 16|, UM, [l8|, which range between 300 and 3000 J/m 2 . The decrease of 
the average shear strength with decreasing fracture energy is explained by the damage 
distribution in the ITZ along the fibre, which is presented together with the shear stress 
distribution in Fig. [9l 



CONCLUSIONS 
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Figure 8: Normalised average shear stress versus embedded length for mechanical analysis 
with the Rankine criterion (0 = 90°) for fracture energies Gf = 30, 300 and 3000 J/m 2 . 
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Figure 9: Contour plot of (a) damage parameter uj and (b) shear stresses in the ITZ 
at maximum pullout force for a frictional angle of (f) = 90° and an embedded length of 
I = 100 /im for a mechanical analysis without the effect of cooling down. 
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In the present work, the influence of radial stresses and the fracture process on the average 
shear strength determined by the microbond test were investigated by nonlinear finite 
element analysis. The work resulted in the following conclusions: 



• Radial stresses generated by the knife blades and the droplet geometry have a small 
influence on the average shear strength. 

• Radial stresses, due to the initial cooling down of the matrix material, increased the 
average shear stress by up to 25 %. 

• Fracture energy influences the average shear stress only for values which are signif- 
icantly smaller than measured in experiments. 



The modelling approach used in this work is based on many simplifications. Material 
nonlinearities are limited to the ITZ between fibre and matrix. All materials are described 
by linear elasticity assuming small strains. Furthermore, the pressure due to the knife 
blade is described by a kinematic constraint at a single node. Nevertheless, the present 
study might help to design future experimental programmes and could provide useful 
guidance for the development of modelling approaches for fibre reinforced composites. 
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